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Abstract 

This is an expanded version of a lecture given at the Workshop on Theoretical Methods for 
Strongly Correlated Fermions, held at the Centre de Recherches Mathematiques, in Montreal, 
from May 26 to May 30, 1999. After general comments on the relevance of field theory to con- 
densed matter systems, the continuum description of interacting electrons in ID is summarized. 
The bosonization procedure is then introduced heuristically, but the precise quantum equiva- 
lence between fermion and boson is also presented. Then the exact solution of the Tomonaga- 
Luttinger model is carried out. Two other applications of bosonization are then sketched. We 
end with a quick introduction to non-Abelian bosonization. 
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1 Quantum Field Theory in Condensed Matter 

Before embarking upon a technical description of what bosonization is and how it helps understand- 
ing the behavior of interacting electrons in one dimension, let us make some general comments on 
the usefulness of quantum field theory in the study of condensed matter systems. First of all, what is 
a quantum field theory? From the late 1920's to the 1980's, it has taken over half a century before a 
satisfactory answer to this fundamental question creep into the general folklore of theoretical physics, 
the most significant progress being the introduction of the renormalization group (RG). A field the- 
ory is by definition a physical model defined on the continuum as opposed to a lattice. However, 
quantum fluctuations being at work on all length scales, a purely continuum theory makes no sense, 
and a momentum cutoff A (or another physically equivalent procedure) has to be introduced in order 
to enable meaningful calculations. Such a regularization of the field theory invariably introduces a 
length scale A~^, which is part of the definition of the theory as one of its parameters, along with 
various coupling constants, masses, and so on. A change in the cutoff A (through a trace over the 
high-momentum degrees of freedom) is accompanied by a modification of all other parameters of 
the theory. A field theory is then characterized not by a set of fixed parameter values, but by a RG 
trajectory in parameter space, which traces the changing parameters of the theory as the cutoff is 
lowered. 

Most important is the concept of fixed point, i.e., of a theory whose parameters are the same 
whatever the value of the cutoff. Most notorious are free particle theories (bosons or fermions), 
in which degrees of freedom at different momentum scales are decoupled, so that a partial trace 
in a momentum shell has no impact on the remaining degrees of freedom. Theories close (in a 
perturbative sense) to fixed points see their parameters fall into three categories: relevant, irrelevant 
and marginal. Relevant parameters grow algebraically under renormalization, irrelevant parameters 
decrease algebraically, whereas marginal parameters undergo logarithmic variations. It used to be 
that only theories without irrelevant parameters were thought to make sense - and were termed 
renormalizable - for the following reason: If the momentum cutoff is taken to infinity, i.e., if the 
starting point of the renormalization procedure is at arbitrarily high energy, then an arbitrary 
number of irrelevant parameters can be added to the theory without measurable effect on the low- 
energy properties determined from experiments. Thus, the theory has no predictive power on its 
irrelevant parameters and the latter were excluded for this reason. In a more modern view, the 
cutoff A does not have to be taken to infinity, but has some natural value Aq, determined by a 
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more microscopic theory (maybe even another field theory) which eventually supersedes the field 
theory considered at length scales smaller than A~^. In the physics of fundamental interactions, 
the standard model may then be considered an effective field theory, superseded by something like 
string theory at the Planck scale (~ 10^'^^cni). Likewise, effective theories for the strong interactions 
are superseded by QCD at the QCD length scale (~ 10~^^cm). In condensed matter physics, the 
natural cutoff is the lattice spacing (~ 10~*cm). 

In practice, field theories should not be pushed too close to their natural cutoff Aq. It is expected 
that a large number (if not an infinity) of irrelevant couplings of order unity exist at that scale, and 
the theory then loses all predictive power. The general practice is to ignore irrelevant couplings 
altogether, and this is credible only well below the natural cutoff.|^ The price to pay for this 
reduction in parameters is that the finite number of marginal or relevant parameters remaining 
cannot be quantitatively determined from the underlying microscopic theory (i.e. the lattice model). 
However, the predictions of the field theory can (in principle) be compared with experiments and 
the parameters of the theory be inferred. After all, this is what the standard model of elementary 
particles is about. But even in cases where a quantitative determination of all the parameters of 
the theory is not possible, some universal (i.e., unaffected by the details of the microscopic theory) 
predictions are possible, and are the main target of fields theory of condensed matter systems. 

The connection between condensed matter systems and relativistic field theories is particularly 
fruitful in one dimension, in good part because the finite electron density may be swept under the 
rug (or rather, the Fermi sea), leaving low-energy excitations enjoying a Lorentz-like invariance in 
the absence of interactions (relativistic field theories are generally applied at zero density). This 
is not true in higher dimensions, where the shape of the Fermi surface makes the connection with 
relativistic field theories more difficult. The development of bosonization was pursued in parallel 
by condensed matter and particle physicists. The former had the one-dimensional electron gas in 
mind, whereas the latter were first interested in a low-dimensional toy model for strong interactions, 
and later by the applications of bosonization to string theory. Ref. reproduces many important 
papers in the development of bosonization, along with a summary of the method (in a language 
more appropriate for particle theorists) and a few historical remarks. 

2 A word on conformal symmetry 
2.1 Scale and conformal invariance 

Let us briefly consider the action of the renormalization group (RG) on fcrmions. The low-energy 
limit is defined around the Fermi surface. The RG action consists in progressively tracing over 
degrees of freedom far from the Fermi surface, towards the Fermi surface. In one dimension, this 
process may be understood better by considering Fig. ^ below. After an RG step wherein the degrees 
of freedom located in the momentum shell [A — dA, A] have been traced over, one may choose to 
perform an infinitesimal momentum scale transformation fc ^ (1 -|- dA/A)fc which brings the cutoff 
back to its initial value A before this RG step. This rescaling is necessary if one wants to compare 
the new coupling constants with the old ones (i.e. if one wants to define a RG trajectory) since we 
are then comparing theories with identical cutoffs A. In a (nearly) Lorentz invariant theory, it is 
natural to perform this scale transformation both in momentum and energy (or space and time). 

In the context of one-dimensional electron systems, Lorentz invariance is an emerging symmetry, 
valid in the low-energy limit, when the linear approximation to the dispersion relation is acceptable. 
The Fermi velocity v then plays the role of the velocity of light: it is the only characteristic velocity 
of the system. Interactions may violate Lorentz invariance and cause the appearance of two char- 
acteristic velocities Vg and Vc (spin-charge separation). However, again in the low-energy limit, the 
system may then separate into disjoint sectors (spin and charge), each benefiting from Lorentz invari- 

^In some applications, leading irrelevant parameters are physically important and must be kept. 
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ance, albeit with different "velocities of light" . In the imaginary-time formalism, Lorentz invariance 
becomes ordinary spatial-rotation invariance. 

A field theory with rotation, translation and scale invariance is said to have conformal invariance. 
Fixed point theories, such as the free boson or the free fermion theories (and many interacting 
theories) are thus conformal field theories. A decent introduction to conformal field theory is beyond 
the scope of such a short tutorial. The subject is quite vast and has ramifications not only in the 
field of strongly correlated electrons, but also in the statistical physics of two-dimensional systems, in 
string theory and in mathematical physics in general. We will simply mention minimal implications 
of conformal symmetry in bosonization. Of course, departures from fixed points break this symmetry 
partially, but their effect at weak coupling may be studied in the conformal symmetry framework 
with profit. 



2.2 Conformal transformations 

By definition, a conformal transformation is a mapping of space-time unto itself that is locally 
equivalent to a rotation and a dilation. In two dimensions, and in terms of the complex coordinate 
z, it can be shown that the only such transformations have the form 

z ad — bc= 1 (1) 

cz + d ^ ' 

where the parameters a, 6, c, d are complex numbers. In dimension d, the number of independent 
parameters of such transformations is {d + l)(d + 2)/2, and d = 2 is no exception. However, the 
case of two dimensions is very special, because of the possibility of local conformal transformations, 
which have all the characters of proper conformal transformations, except that they are not one- 
to-one (they do not map the whole complex plane onto itself). Any analytic mapping z ^ w on 
the complex plane is locally conformal, as we known from elementary complex analysis. Indeed, the 
local line element on the complex plane transforms as 

The modulus of the derivative embodies a local dilation, and its phase a local rotation. This 
distinctive feature of conformal symmetry in two-dimensional space-time is what allows a complete 
solution of conformal field theories, even in circumstances that apparently break scale invariance. For 
instance, the entire complex plane (the space-time used with imaginary time at zero temperature) 
may be mapped onto a cylinder of circumference L via the complex mapping z = ^'^'^1^ , This 
allows for the calculation of correlation functions in a system with a macroscopic length scale (a 
finite size L at zero temperature, or a finite-temperature /3 = L/w in an infinite system) from the 
known solution in a scale-invariant situation. Mappings may also be performed from the upper 
half-plane unto finite, open regions; such mappings are particularly useful when studying boundary 
or impurity problems. On the formal side, this feature of conformal field theory makes the use of 
complex coordinates very convenient: In terms of x and imaginary time r = it, these coordinates 
may be defined as 

z = —i{x — vt) = VT — ix , . 

z = i{x + vt) = VT + ix ^ ' 

with the following correspondence of derivatives: 

d-z = -^(^5* + ^-) dt = iv{dz + dz) (4) 
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A local operator 0{z, z) belonging to a conformal field theory is called primary (more precisely, 
quasi-primary) if it scales well under a conformal transformation, i.e., if the transformation 

0'{az, az) = a-^a-~^0{z, z) (5) 

is a symmetry of correlation functions. The constants h and h are a fundamental property of the 
operator O and arc called the right- and left-conformal dimensions, respectively. Under a plain 
dilation (a = a), the operator scales as 0'{az, az) = a~^0{z^z), where A = ft, + ft, is the ordinary 
scaling dimension of the operator. Under a rotation (or Lorcntz transformation), for which a = e'^ 
and a = c~*^, the operator transforms as 0'{e^^z, e~^^z) = e^^^O{z,z), where s = ft — ft is called 
the conformal spin of the operator. 

The scaling dimension appears directly in the two-point correlation function of the operator, 
which is fixed by scaling arguments - up to multiplicative constant: 

iO{z,z)OH0,0)) = ^S- (6) 



2.3 Effect of perturbations 

Fixed-point physics helps understanding what happens in the vicinity of the fixed point, in a per- 
turbative sense. Consider for instance the perturbed action 

S = So + g j dtdx 0{x, t) (7) 

where is the fixed-point action and O an operator of scaling dimension A. Let us perform 
an infinitesimal RG step, as described in the first paragraph of this section, with a scale factor 

A = 1 -t- dA/A acting on wavcvcctors; the inverse scaling factor acts on the coordinates, and if x 
and t now represent the new, rescaled coordinates, with the same short-distance cutoff as before the 
partial trace, then 

S' = So + g j d{Xt)d{Xx) 0{Xx, Xt) 

= So + gX^-^ j dtdx 0{x, t) (8) 

We assumed that we are sufficiently close to the fixed point so that the effect of the RG trace in 
the shell dA is precisely the scaling of O with exponent A. By integrating such infinitesimal scale 
transformations, we find that the coupling constant g scales as 



/ . \ A-2 / /■ \ 2-A 



90 V C 



(9) 



where ^ and arc correlation lengths at the two cutoffs A and Aq. Thus, the perturbation is relevant 
if A < 2, irrelevant if A > 2 and marginal if A = 2. 

A relevant perturbation is typically the source of a gap in the low-energy spectrum. In a Lorentz- 
invariaut system, a finite correlation length is associated with a mass gap m ^ Let us suppose 

that the correlation length is of the order of the lattice spacing ~ 1) when the coupling constant g 
is of order unity. The above scaling relation then relates the mass m to the bare (i.e. 'microscopic') 
coupling constant: 

l/(2-A) /T„x 

mr^go" (10) 

Given the warning issued in the introduction on the impossibility of knowing the bare parameters of 
the theory with accuracy, we should add that 'bare' may mean here "at a not-so-high energy scale 
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where the coupUng go is known by other means". In fact, this scahng formula works surprisingly 
well even if go is taken as the bare coupling (at the natural cutoff Aq), provided go is not too large. 

Conformal field theory techniques may also be used to treat marginal perturbations, in perform- 
ing the equivalent of one-loop calculations in conventional perturbation theory, with the so-called 
operator product expansion. This is explained in Appendix A and an example calculation is given 
(more can be found in Refs ||^, || ^). However, perturbations with conformal spin s ^ are more 
subtle, and typically lead to a change in kp ot incommensurabilities. 

2.4 The central charge 

A conformal field theory is characterized by a number c called the conformal charge. This number 
is roughly a measure of the number of degrees of freedom of the model considered. By convention, 
the free boson theory has conformal charge c = 1. So does the free complex fermion. A set of N 
independent free bosons has central charge c = N . Conformal field theories with conformal charge 
c < 1 correspond to known critical statistical models, like the Ising model (known since Onsager's 
solution to be equivalent to a free Majorana fermion), the Potts model, etc. Models at c = 1 
practically all correspond to a free boson or to slight modifications thereof called orbifolds, and have 
been completely classified . A conformal field theory of integer central charge c = N may have a 
representation in terms of N free bosons, although such representations are not known in general. 
Likewise, a conformal field theory of half-integer central charge may have a representation in terms 
of an integer number of real (i.e. Majorana) fermions. 

3 Interacting electrons in one dimension 

3.1 Continuum fields and densities 

Let us consider noninteracting electrons on a lattice and define the corresponding (low-energy) field 
theory. This poses no problem, since the different momentum scales of this free theory are decoupled. 
The microscopic Hamiltonian takes the form 

Hf = ^e{k)c''{k)c{k) (11) 

k 

where c(k) is the electron annihilation operator at wavevector k (we ignore spin for the time being). 
The low-energy theory is defined in terms of creation and annihilation operators in the vicinity of 
the Fermi points (cf. Fig. |l|) as follows: 

a{k) — c{kp + k) 

a{^k) = c{—kp — k) 

Pik) ^ c\kp-k) 

I3{~k) = c){-kp + k) (12) 

where k is positive. The noninteracting Hamiltonian then takes the form 

ffp = y '^v\k\ {aHk)a{k) + p\k)m) (13) 

where v is the Fermi velocity, the only remaining parameter from the microscopic theory. The energy 
is now defined with respect to the ground state and the momentum integration is carried between 
—A and A. The operators a{k) and /3(fe) respectively annihilate electrons and holes around the right 
Fermi point {k > 0) and the left Fermi point (A: < 0). A momentum cutoff A is implied, which 
corresponds to an energy cutoff vA. Note that the continuum limit is in fact a low-energy limit: 
Energy considerations determine around which wavevector (±kp) one should expand. In position 



6 




Figure 1: Typical tight-binding dispersion in ID, illustrating left and right Fermi points and the 
linear dispersion in the vicinity of those points. 



space, this procedure amounts to introducing slow fields and such that the annihilation operator 
at site n isQ 

= V(a;)e''^'^^+V^(a;)e-*'=^^ (14) 

The factor ^/a is there to give the fields the proper delta-function anticommutator, and reflects 
their (engineering) dimension: 

{V;(a;),i/;^(a;')} = S{x - x') 
yix),^pf{x')} = S{x-x') 

{^ix),i^\x')} = (15) 



Left-right separation 

The mode expansions of the continuum fields are 

Jk>0 

Hx) = I '^[e^^^a{k)+ e-^^-p^k)] (16) 

The time dependence of a(fc) and /3(fc) is obtained through multiplying by the phase ^ in 

terms of the complex coordinates (^ , the mode expansions for the time-dependent fields are then 

^{z) = / [e-'=^a(fc)+e'=^/3t(fc)] 

m = I ^[e^^"a(fc)-Ke-'=^"/3t(fc)] (17) 

Thus the right-moving field V' depends solely on the right- moving coordinate z, whereas the left- 
moving field t\} depends solely on z. 

The mode expansion (|l^) is misleading in one respect: it makes believe that the positive- 
wavevector modes are contiguous to the negative- wavevector modes in fc-space, which is not the case 

^We will generally use a bar ( ) to denote left-moving operators, and the same symbol without the bar for right- 
moving operators. A more common notation in condensed matter physics is the use of subscripts L and R. The 
(lighter) notation used here stresses the analogy with complex coordinates. 



7 



(they are separated by roughly 2kp). Thus right- and left-moving modes are clearly well-separated. 
Left and right electrons have density fluctuations 

J{x) = ^^x)^P{x) J{x) = 'ip'iix)ip{x) (18) 

The total electronic density (with respect to the ground state) being rttot. — J + J- 

Analogy with the Dirac equation 

It is a simple matter to express the continuum Hamiltonian (um in terms of the fields ip and ip: 



Hp = -iv / dx [^''d^^ - i^'^d^ij] (19) 



This is exactly the Dirac Hamiltonian. Indeed, in (1+1) dimensions, the Dirac matrices can be taken 
to be two-dimensional, for instance 

/^(° l) l) (20) 



^-1 oy ' V 1 

and the two-component Dirac spinor can be taken to be 



*=(S) (21) 



The Dirac Lagrangian density - with a mass term - being 

Cd = i*(7°5t - W7^ax)* - m-^-^ (22) 
The canonical Hamiltonian density derived from this Lagrangian is therefore 

Hd = -iv [t/j^d^ip - li^^d^Tp] + m{ilj^ - V^^V) (23) 

which coincides with Eq. ( p^ in the massless case. 

Adding a mass term to the continuum Hamiltonian amounts to opening a gap 2m at the Fermi 
level. On the other hand, adding a term like 

n' = ii)H + i>U^ = J + J (24) 

to the Hamiltonian would simply shift the chemical potential, i.e., modify the value of kp. This 
would of course require a redefinition of left and right modes, but the end result would again be a 
Hamiltonian of the form (p^, with a slightly different velocity v. 

In the language of Conformal Field Theory, the free fermion theory has central charge c = 1, and 
the conformal dimensions of the fermion fields are 

: {h~h)^{\,Q) ^ : {h~h) = {Q,\) (25) 

Thus, the mass term ^'^'tp — ip'^^p has scaling dimension A = 1 and zero conformal spin; according 
to Eq. (0), this term gives rise to a gap proportional to m, as expected of course. On the other 
hand, the perturbation tp + : even if it also has scaling dimension 1, is a sum of terms 
with nonzero conformal spin. For such terms the analysis leading to Eq. ( |l0[ ) cannot be applied. 
Perturbations with conformal spin, such as this one, will typically change the slow wavevectors, 
cause incommensurability, etc. 
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Figure 2: The four scattering processes for right-moving (continuous lines) and left-moving (dashed 
lines) electrons in one dimension: (i) 514 is the amplitude for forward scattering; (ii) 172 for dispersion 
(scattering of left onto right electrons); (iii) gi for backward scattering (each electron changes branch, 
but keeps its spin); (iv) (73 for Umklapp scattering (two left electrons become right electrons or vice- 
versa; possible only at half- filling because of crystal momentum conservation). Spin indices are 
suppressed, but incoming electrons must have opposite spins for gi and (73. 



3.2 Interactions 

The general two-body interaction between ID electrons (of one band) is 

-ffint. = J2 V„,„,{kuk2,h)ci{k^)cl,{k2W{k^)c^{ki) (26) 



where we have restored spin indices. In the low-energy limit, scattering processes are naturally 
restricted to the vicinity of the Fermi points, and fall into four kinematic types, illustrated on Fig. ^j. 
In terms of continuum fields, the corresponding Hamiltonian densities are: 



Hi 


= vgi X! V'iV'.rV'-.TV'-a 




-HI 








= vg2,s{JT - Ji){Jt - Ji) 




^3 




-H.c. 


HI 


CJ 

= ^vgi,a[{J^ +Jlf + {J^ 




HI 


= ^vg4,s [iJT - Jif + [J^ 





(27) 
Remarks: 

1. The incoming spins in Hi can be assumed to be antiparallel, since the case of parallel spins 
reduces to a density-density coupling of the form Tij- 

2. We have supposed that the couplings 51-4 are momentum independent, and this can be safely 
assumed in the low-energy limit. However, they are cutoff dependent, i.e., they are subject to 
a RG flow. 

3. The interactions Hi and are not separately spin-rotation invariant, but their combination 



is if gi — ~2g2.s- This will be shown in subsection 7.2 



That 52,/^ and 34, ^ are density-density interactions is made more explicit in this real-space 
representation. If these couplings alone are nonzero (i.e., gi — gs — 0) , the resulting model is 
called the Tomonaga-Luttinger model and can be solved exactly by bosonization. 



The bosonized form of these interactions will be derived in Sect. 5.6. 
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Figure 3: Particle-hole excitations in one dimension. 




Figure 4: Particle-hole excitations in two dimensions. 

4 Bosonization: A Heuristic View 

4.1 Why is one-dimension special? 

The basic idea behind bosonization is that particle- hole excitations are bosonic in character, and 
that somewhat the greatest part (if not the totality) of the electron gas spectrum might be exhausted 
by these excitations. The question was raised by F. Bloch back in 1934, but Tomonaga realized in 
1950 that it could only be true in one dimension. The reason is simple: Consider Fig. ^. On the 
left the creation of a particle-hole pair of momentum k is illustrated. On the right is the spectrum 
of particle-hole excitations created out of the Fermi sea: the energy of the pair is plotted vs the 
momentum of the pair (both quantities measured w.r.t. the ground state). Because of the linear 
one-particle dispersion near the Fermi level, the pairs have a narrow, quasiparticle-like dispersion 
near zero momentum: they can propagate coherently. In other words, the particle and hole have 
nearly the same group velocity at low energies and propagate together. Any weak particle-hole 
attraction is then bound to have dramatic effects, i.e., to bind the pair into a coherently propagating 
entity: a new particle. 

On Fig. ^, the corresponding situation in two dimensions is illustrated, with a circular Fermi 
surface. There it is clear that a particle-hole pair with a given momentum k can have a continuous 
spectrum of energies, starting from zero. Thus, the particle-hole spectrum is a continuum throughout 
and interactions have a harder time forming coherently propagating particle- hole pairs In any case, 
defining the theory in terms of bosonic excitations is much less obvious in dimensions greater than 
one. 

In the low-energy limit, particle- hole excitations in on dimension are exactly coherent and thus the 

^This argument is partially weakened by the presence of nesting. 
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density fluctuation fields and J„ create propagating particles. The question is then whether any 
operator can be expressed in terms of J^- and J„ , especially a single-particle creation or annihilation 
operator. The answer is yes; the right-moving fermion field may actually be expressed as 



1 



'2tt 



exp 



2m 



Ax' J Ax') 



(28) 



The rest of this section, and much of the following one, will be devoted to proving Eq. (|7^), equivalent 
to the above. 



4.2 The simple boson 

Before going any further one must define what a free boson is. This may be done fundamentally 
from a Lagrangian. For a massless boson field (f, this is 



^0 = ^ / 



V 



(29) 



where v is the speed of the Bose particles. The corresponding Hamiltonian density is then 

Ho - Iv [tf + {d^ipf] n = -dtip (30) 
where 11 is the field conjugate to ip, the two fields obeying the canonical commutation rules 

Mx)Mxl]^^Six-x') -J (31) 

The standard mode expansion for the fields ip and 11 on the infinite line is 

dfc I V r. ,kx , Lt^;,\ „-ifcxl 



27r V 2uj{k) 



e 



n(x) = J^^^[-,b{k)e^'^+^bHk)e-^'^^] (32) 

where uj{k) = v\k\. The creation and annihilation operators obey the commutation rules 

[b{k),b\k')]^2n6{k^k') (33) 
and the Hamiltonian may be expressed as 

/dk 
—Lo{k)b\k)b{k) (34) 

The ground state (or vacuum) |0) is annihilated by all the b{k)'s, and has zero energy by convention. 

Like fermions, Bose particles in ID cither move to the left or to the right. Indeed, the time 
dependence of the annihilation operators is obtained through multiplying by the phase e"*"''''*, if 
the boson is massless. A separation into left- and right-moving parts is then possible. We write 

if{x,t) ^ (f){x ~ vt) + (j){x + vt) (35) 

where 

0(z)= / [6(fc)e-'=^+6t(fc)e'=^] 



/fe>o 2ny/2k 

m = / ^4= \b{k) e-^^' + b\k) e'^^l (36) 
Jfe>o '^^ V 2k 

The variables z and z are defined as in Eq. (||). Note that wavevectors now take only positive values, 
since we have defined b{k) = b{—k). 
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The dual field 



It is customary to introduce the so-called dual boson d^x^t), defined by the relation dxi^ = —11 = 
— ^dtf- In terms of right and left bosons (p and 0, this becomes (cf Eq. (^) 

V 

= -i{dz - ds)i(f> - <j)) 

= dx{cf>-^) (37) 
therefore 9 — (j) — (j), modulo an additive constant which we set to zero. One may then write 

0-1(^ + 7?) 0^i(^-,9) (38) 

if one desires a definition of (j) and (j) that stands independent from the mode expansion. Since "d is 
expressed nonlocally in terms of 11 (i.e., through a spatial integral), the above expression makes it 
clear that the right and left parts 4> and cj) are not true local fields relative to (p. The basic definition 
of •& in terms of 11 implies the equal-time commutation rule 

[ip{x),d{x')]== ~i9{x-x') (39) 

(please note that 9{x) is the step function [7^ d]) which further shows the nonlocal relation between 
if and ^ (the commutator is nonzero even at large distances). 



4.3 Bose representation of the fermion field 



After this introduction of the boson field, let us proceed with a heuristic (and incomplete) derivation 
of the bosonization formula. The basic idea is the following (we ignore spin for the moment): the 
electron density ntot. being bilinear in the electron fields, it has Bose statistics. Let us suppose that 
it is the derivative of a boson field: ntot. oc d^^- Then 



Lp{x) = A / Ay ntot.iy) 



"■tot. (a;) = -jdxfix) 



(40) 



where A is a constant to be determined. Creating a fermion at the position x' increases by A if 
x < x' . This has the same effect as the operator 



exp 



-iX 



dy n(y) 



(41) 



which acts as a shift operator for ip{x) if x < a;', because of the commutation relation 



ip{x), 



dy n{y) 



{x < x') 
{x > x') 



(42) 



Thus, a creation operator ^p^{x') or ^p^{x') may be represented by the above operator, times an 
operator that commutes with (p {e^°"^, for instance). This extra factor must be chosen such that -0^ 
depends on x — vt only, and on x + vt only. The left-right decomposition ( ^ ) and Eq. ( ^ ) are 
useful here. Since 

n = -dtp) = -dx^ (43) 
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we have 



r' 

I dy n{y) = iXiiix') (44) 



from which purely right- and left-moving fields may be obtained by adding or subtracting —^iaLp{x'). 
Thus, a sensible Ansatz for a boson representation of electron creation operators is 

V'^ (x) = A e2*^'^(^) V'^ [x) = A e'-^i^^i^) (45) 

where the constants A and A are to be determined, for instance by imposing the anticommutation 
relations (|l5|). This will be done below. 

The one-electron states V'^lO): being generated by exponentials of boson creation operators, are 
coherent states of the boson field (j)^ On the other hand, the elementary Bose excitations (the mesons, 
in field-theory jargon) are collective density (or spin) fluctuations, i.e., electron-hole excitations. 

5 Details of the Bosonization Procedure 
5.1 Left and right boson modes 

The mode expansions (^) may seem reasonable, but are actually incomplete, for the following 
reasons: 

1. The field being massless, the mode expansion is ill-defined at fc = 0, because uo{k = {)) = 0. 
The mode at fc = 0, or zero-mode, has to be treated separately and turns out to be important. 
The existence of this mode apparently spoils the left-right separation of the field Lp. 

2. In order for bosonization to be rigorously defined, the field must have an angular character. 
In other words, the target space of the field must be a circle of radius R (to be kept general 
for the moment), so that Lp and Lp + 271 R are identified. We say that the boson is compactified 
on a circle. 

3. A rigorous proof of bosonization procedure (at the spectrum level) is best obtained on a system 
of finite size L with periodic boundary conditions. 

Taking the above remarks into account leads to the following improved mode expansion for the 
boson field Lp{x, t): 

^{x, t)=q+^ + ^ + J2^[bn e''^ + bi e^^- + K e-^^- + H e^'^] (46) 



n>0 ^ 

Wavevectors are now quantized as fc = 21:71/ L (n an integer) . The creation and annihilation operators 
have a different normalization from Eq. (|32|), so as to obey the commutation rules 

[&„, 6|„] = (5m„ [&n,6j,J = (5m„ (47) 

The zero- mode is treated explicitly as a pair of canonical variables q and ttq obeying the commutation 
rule 

[q, TTo] = i (48) 
Defining 1^9 on a circle of radius R has two consequences: 

^The interactions gi and 93 give the boson Hamiltonian the sine-Gordon form (cf. Sect. 5.6). In that case, 
one-electron states correspond to soUtons (kinks) of the sine-Gordon theory. 
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1. The operator q is not single- valued, since it is an angular variable {q and q+2T:R are identified). 
Only exponentials e*"'?/-'^, where n is an integer, are well defined. Because ttq is conjugate to 
q, we have the relation 

^-^nq/R^^ ^inq/R ^ 7^0 + - (49) 

R 

(the exponential operators shifts the eigenvalues of ttq). Starting with the ttq = ground state, 
this means that the spectrum of ttq is restricted to Z/i?. 

2. Winding configurations are allowed; the constant ttq is precisely defined as ttq = 2T:Rm, where 
TO is the number of windings oiip as x goes from to L. It turns out to be useful to introduce an 
operator ttq, defined by its eigenvalues 27r i?m. This operator commutes with all other operators 
met so far, but one may define an operator q conjugate to ttq, i.e., such that [g, ttq] = i. 

The separation between left- and right-moving parts may then be done in spite of the existence 
of a zero mode: 

P , , ^ 1 , 



Q + :77(wt - a;) + V -j= (&„ 6"'=" -t- 1 

n>0 ^ 



(z) = Q+f^{^vt + x) + Y,^=(Ke-'^' + hie^') (50) 



2i -^-^ v47m 

where we have defined the right and left zero-modes 

Q = \{<1-S[) P = 7:o-n„ [Q,P] = i 

Q = ^(q + q) P-^o+^o [Q,P]^i (51) 

Note that the "artificial" operator q drops from the combination ip = (f> + (j). The spectrum of P and 
P is 

P= + 2TTRm^ P = - 27rPTO) (52) 

TO and n being integers. Note that this spectrum is not a simple Cartesian product, which means 
that the left-right decomposition of the boson is not perfect, despite the expansion (pC^ 



5.2 Proof of the bosonization formulas: Vertex operators 

Before investigating the constants A and A of Eq. (^) , serious questions must be asked on the precise 
definition of an exponential operator of the form e'"'^ (such operators are called vertex operators, 
from their usage in string theory). Because (f) is a fluctuating field, any power or more complicated 
function of (j) requires a careful definition if it is not to have a divergent average value. The natural 
prescription is called normal ordering and consists in expressing (f> as a mode expansion and putting 
all annihilation operators to the right of creation operators. Thus, the exponential operator has the 
following precise definition: 



ia bj e'^^ exp ia 



n>0 



n>0 



-iaP(x-vt)/2L 



(53) 



where P is considered an annihilation operator since the boson ground state |0) has zero momentum. 
The notation :A: means that the operator A is normal ordered. 

Because of normal ordering, vertex operators do not multiply like ordinary exponentials. Instead, 
one has the relation 

(54) 
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where {4i{z)4i{z')) is a simple ground state expectation value if the left-hand side is an ordinary 
product, or a Green function if it is a time-ordered product. This is one of the most important 
formulas of this review, and it can be easily demonstrated using the Campbell-Baker-Hausdorff 
(CBH) formula: 

e^ = e^+^ e^I^'^1 ([A, B] = const.) (55) 
Consider for instance a single harmonic osciUator a, and two operators 

A = aa + a' a) B = f3a + (56) 

The CBH formula allows for a combination of the normal-ordered exponentials: 

■ ■ ■ e^ ■ = e"'"* e"" e^'"'^ e^°' 

= e"'"* e'^'"* e"" e'^" e"'^' 

The last equality also applies if (a, a^) is replaced, in the same order, by a pair (p, q) of canonically 
conjugate operators ([(7,p] = «)• It also applies to a combination of independent oscillators, such 
as the boson fields <j) or (j). Finally, it applies to a time-ordered product as well as to an ordinary 
product, as long as (0|A_B|0) is also replaced by a time-ordered product. Thus, the identity (|^ is 
proven. 

It remains to calculate the Green function {(p{z)(p{z')). This can be done in several ways. For 
instance, by using the mode expansion (|50|) and taking the limit L — > oo: 



We will adopt the normalization 



2L ' ' An ^ n 

n>0 

Z 1 ^ 2'KZ 

In z -I- const. (L — > cx)) (58) 

A-K 



{cp{z)^{z'))^-^Hz-z') (59) 



Att 

and drop the constant term, which basically defines an overall length scale. Likewise, one finds 

{^{-z)^{-z')) = -l-\n{-z--z') (60) 

Another way of computing the boson Green function is to notice that the nonchiral Green function 
Q{x^t) = ((/3(a;, r)(y9(0, 0)) must obey the two-dimensional Poisson equation 

V'^g = 5{x)5{vt) (61) 

in Euclidian space-time. The solution to that equation is readily obtained in polar coordinates, up 
to an additive constant (the correct normalization is obtained by applying Gauss' theorem around 
the origin): 

G{x,vt) = -^\n{x^ +v^t'^) = --^-ln(zz) (62) 
At: Att 

The result is thus equivalent to what can be inferred from the mode expansion, since 

(^(x, tMO, 0)) - (^(z)0(O)) + (mm) (63) 
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Formula ( p^ can thus be rewritten as 

We are now in a position to demonstrate the boson-fcrmion equivalence at the level of the 
(anti)commutation relations. In fact, it is simpler to determine the constants A and A of Eq. ( |45| ) by 
comparing Green functions rather than simply looking at equal time (anti)commutators (this avoids 
the equal time singularity). The electron propagator is readily calculated from the mode expansion 



fc>0 Jq>0 27r 
^Q-k{z-z') 

fe>o 27r 

^ (65) 



2tt z — z 

where we have supposed that t > t', which garantees the convergence of the integral and follows 
from the prescribed time-ordered product. Likewise, 

(V.t(z)^(z')) = ^^ (66) 
The correct anticommutator may be obtained from the equal time limit: 

{i;{x,0),ip\x',0)} = \im({iP{e-ix)^p''{-ix')) + {ij^{e-ix')ip{-ix))) 

£—►0 

lim 



27r e^o ye — i[x — x') e — i{x' — x) 
= ^^^^7 I 2 = H^-x') (67) 

(■0 or tp appearing with a single argument are considered functions of z or z). One easily checks that 
the space-time Green function (Bq) has the expected expression in momentum-frequency space: 



1 1 



(68) 



2TT z — z' uj — v\k\ 
To demonstrate Eq. (^) and fix the values of A and A, one needs to calculate 

{^{z)ij\z')) = A2(e~2«A0(^) ^2^X,piz')^ ^gg) 

From the relation (|6^), this is 

{iPiz)ij^iz')) = A2(e-2'M.^'W-^("'»)(z-z')"^'/" 

= A^z- z')-^'/'' (70) 

Note that the normal-ordered expectation value 

f^^-ia(t>{z)-il3(t>(z')^ (-JX) 

is nonzero only ii a + /3 = 0. This can be inferred from the action of the zero- mode operator Q:^ 

(0|e»("+««e-*("+'^)^/'^|0) (72) 

^For a calculation in the path-integral formalism, see Ref. j6|. 

^It is also a consequence of the U(l) symmetry ip ^ ip + a, where a is a constant. 
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vanishes if a + /? 7^ 0, since the exponential of Q acting on |0) yields another eigenstate of P, with 
eigenvalue a + f3, orthogonal to |0). The condition a + /? = is called the neutrality condition]^ By 
comparing Eq. ( [70| ) with Eq. (|65|), we conclude that the constant A must be ^/^^ and that A must 
be l/\/27r. The correct bosonization formula is then 

(73) 

Mx) ^^x) =^Le-*^"?(^) 



Notice the change of sign in the exponent between -0 and t/i. Notice also that the choice of additive 
constant leading to Eq. ( |60| ) influences the value of A only, the end result being the same irrespective 
of that choice. 
Remarks: 

1. It is also instructive to check the anticommutation relation 

{V'(x),V^(a;')} = (74) 
This is done simply by noting that the Green function 

{iP{z)xl;(z')) = _L(e-*^^(^)e-*^^(^')) (75) 

vanishes, by the neutrality condition. The anticommutator, obtained in the equal-time limit, 
also vanishes. 

2. The vertex operator e'"*^ is a scaling field. Its conformal dimension may be retrieved from 
the exponent of the Green function 

^^.a0(.)(^,o0(o))t) = _J_ (76) 

and from the relation (|^) . The same analysis may be applied to a left-moving vertex operator 
e^"^. Thus, 

h(a) = ^ h{a) =0 

^ ' ^ ^ ' _2 (77) 

h{a) =0 h{a) =^ 

Roughly speaking, the boson field has scaling dimension zero and its powers series may produce 
operators with nonzero scaling dimension because of the necessary normal ordering. 



5.3 Bosonization of the free-electron Hamiltonian 

The bosonization formulas ( [73| ) are the basis of a fermion-to-boson translation of various operators. 
However, this translation process is often subtle, since normal ordering of these operators is required. 
Take for instance the chiral fcrmion densities J and J. Normal ordering can be constructed from 
the mode expansion. But a shorter and more elegant way to achieve it is by point splitting^ i.e., by 
defining 

J{z) = lini [V/ {z + e)ilj{z) - (V/ {z + e) V'(z))] (78) 

^This terminology comes from an analogy with the Coulomb gas system of two-dimensional statistical mechanics. 
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where e has positive time (so as to appear as shown after time ordering). This hmit must be taken 
after applying the bosonization formulas and the relation (pj): 



J 



— lim 

27r e^o 



— lim 

2-K e^O 



— lim 

27r £^0 



^ieV4^d,4'iz)l, _ 1 

e e 



2tt 



(79) 



where we have Taylor-expanded the exponential on the last line. We proceed likewise for the left- 
moving density, except that the sign is opposite. To summarize: 



J 



J = 



zdzf 



(80) 



(note that 9^0 — dz(p and 9^0 = dztp). 

The density fluctuations J and J are also components of a conserved current from the boson 
point of view. Indeed, the compactified boson has a U{1) symmetry: the Lagrangian is invariant 
under the shift (p — > ip + a. By Noether's theorem, this symmetry implies the existence of a conserved 
quantity, with a density proportionnal to dt^p and current proportionnal to v^d^ip. It is customary 
to define right and left 'currents' exactly as in Eq. (|8^), so that the continuity equation reduces to 
dzJ + dzJ = 0. In fact, the critical (i.e., massless) and one-dimensional nature of the theory enhance 
this U{1) symmetry to a chiral U{1) symmetry, by which the left and right currents are separately 
conserved: 

dzJ = dzJ = (81) 

The same short-distance expansion as in Eq. ( |79| ) may be used to demonstrate explicitly the 
equivalence of the electron Hamiltonian (|l^) with the boson Hamiltonian ( ^0[ ) , except that a higher- 
order expansion is necessary. More explicitly, an expression like ip^dx4' must be evaluated as 



-i hm {^^{z + e)dz^{z) — {ip^ {z + e)dz^{z))] 

e— 



The limit is evaluated by using the following short-distance products (to order e^) 
bosonization formulas and the relation (M): 



(82) 

derived from the 



^{z')i:^{z) = 


1 

27re 


+ i—;=dzip + i- 




e{dz^f 


^Hz'Mz) = 


1 

27re 


I ^dz'P I 




e{dzv? 


^P{z')^Hz) - 


1 

27ri 


I ^dzf I 




e{d-zvf 




1 

27ri 


-1- i^=dzf + i- 




e{d-zvf 



(83) 



where e 



and e = z' — z, and the fields on the r.h.s. are considered functions of (z, z) only. 



the (z', z') dependence residing in the powers of e and e only. The limit ( |82| ) may be obtained by 
differentiating with respect to z or z the above expansions, and proceeding to a straightforward 
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substitution. The Hamiltonian (19) then becomes 



= -V j Ax [{d^ipf + {d,vf] 

= TTv [ dx ( J2 + P) (84) 



which is precisely the boson Hamihonian (|3^), after using the relations (^. Thus the equivalence 
between a free fermion and a free boson is demonstrated at the level of the Hamiltonian. Note that 
we have not demonstrated this equivalence at the Lagrangian level; this cannot be done with the 
techniques used here. 

5.4 Spectral equivalence of boson and fermion 

That the free boson (at a certain radius R) is equivalent to a free fermion requires that the spectra of 
the two theories be identical. This can be tested by comparing the grand partition functions of the 
two models. Let us try to do it for left and right excitations separately. Consider the right-moving 
fermion "tp^x), and the corresponding Hamiltonian and fermion number: 

i/F =^^^A:(4afc+4/3fe) iVp = ^(4"^ " ^/^fc) (85) 

fe>0 fc>0 

where the operators ak and Pk have different normalizations from those of Eq. ([l3| ) and obey the 
anticommutation rules 

{ak, ckg} = Skq {ak,aq} = (86) 
(and likewise for f3k). The spectrum of the Hamiltonian is encoded in the grand partition function 

Zf=Y^ e-'^f^-^^) (87) 

states 

We have defined the theory on a cylinder, so as to have discrete wavenumbers. We shall use an- 
tiperiodic boundary conditions, known in string theory as Neveu-Schwarz boundary conditions: 

iPix + L) = -i^ix) (88) 

The wavenumbers are then half-integer moded: k = 27r(n + \)/L, where n is a positive integer (by 
contrast, periodic boundary conditions, also known as Ramond boundary conditions, lead to integer- 
moded wavenumbers and to the presence of a fermion zero- mode). Each mode being independent, 
the fermion grand partition function factorizes as 

oo 

-Zf = n (1 + + g"+^/^<"^) q = e-2'^'"^/^ , t ee e^'' (89) 

n=l 

Let us now turn to the boson. We first have to write down an expression for the Hamiltonian in 
terms of the modes appearing in (|5^) . A straightforward calculation yields 

^B = — + — ^ (90) 

n>0 

for the right-moving part. Let us then notice that the fermion number N-p is proportional to the 
zero-mode P. More precisely, 

iVp - ^ (91) 
\/47r 
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This can be shown from the mode expansion ( pOD and by the equivalence J = —dxipl \/t^, which 
implies 

P 



Jdxj = -^WL) - m) = 



/An 

Each mode &„ being independent, the boson grand partition function is then 



(92) 



I n (1 + 9" + '?'" + '?'" + ■■•)[ E q^'^^^^tP/^ 

[n=l J P 



(93) 



So far we have not determined the boson radius R, nor the spectrum of P. We now invoke Jacobi's 
triple product formula: 

oo 

[|(l-g")(l+(Z"-i/2^)(l + g"-i/2t-i) = ^9"'/2t" (94) 

n=l neZ 

Applied to the fermion partition function, this remarkable identity shows that 



l,n=l ^ ) n 



(95) 



In order to have Zp = Zb, we therefore need P to take the values P = n^/An. This is not quite the 
spectrum (|5^) . The closest it comes to is obtained at radius R= 1/ V47r, where (Hj) becomes 

P = \/47r(n + im) P = V^{n - ^m) (96) 



In fact, because of the imperfect left-right separation of the boson spectrum (|52|), the equivalence 
with a perfectly left-right separated fermion theory is impossible. However, a state-by-state corre- 
spondence between the two theories exists if some constraints are imposed on the fermion spectrum. 
Explicitly: 

1. One needs to include both periodic and antiperiodic boundary conditions in the fermion theory: 
the Hilbert space is then the direct sum of two sectors. 

2. The periodic sector must contain a single zero- mode ■00 (the same for left- and right- moving 
fermions), such that {tpo, V'q} = 1- The operator ipl does not change the energy of a state, but 
increases the fermion number by one. This mixes (in a weak way) the left and right-moving 
spectra. 

3. the number of fermions in the antiperiodic sector must be even (odd- number states are thrown 
out). 

These constraints will not be demonstrated here; a detailed proof can be found in Refs and [Q. 

That the simple-minded correspondence between left-moving fermions and left-moving bosons 
is impossible should not bother us too much. The complications mentioned above basically involve 
boson or fermion zero-modes and boundary conditions. The nonzero frequency modes are not 
affected and the fcrmion-boson correspondence works well for those modes, as reflected in the Jacobi 
triple product formula (p4). First and foremost, we should retain from this exercise that the boson 
radius must be i? = l/y/An in order for the fermion-boson equivalence to hold. 
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5.5 Case of many fermion species: Klein factors 



Suppose now that we have more than one fermion species, labelled by a greek index: -0^. Bosoniza- 
tion would then require an equal number of boson species and the formulas ( [73| ) would still be 
applicable, except for the fact that they do not provide for the anticommutation of different fermion 
species. In fact, this problem arises for a single species as well, since the left and right fermions 
and '0 must anticommute, whereas the fields 4> ^-nd (j) commute.^ A possible solution lies in the 
introduction of additional anticommuting factors, the so-called Klein factors, in the bosonization 
formulas: 



1 



27r 
1 



1 



2tt 
1 



The Klein factors r/^ and rj^ are Hermitian and obey the Clifford algebra: 



(97) 



(98) 



Klein factors act on a Hilbert space distinct from the boson Hilbert space generated by the modes 
of Eq. (50). This Hilbert space expansion must be compensated by some sort of "gauge fixing": the 
Hamiltonian must be diagonal in this Klein factor space, as well as required physical observables. 
Then, one Klein eigenstate is chosen and the rest of the Klein Hilbert space decouples. 

An alternative to Klein factors, which requires no extra Hilbert space, is to include in the 
bosonization formula for ijj^ the factor 



exp 



= exp 



(99) 



{Pi, is the conjugate to the zero mode Q^, cf. Eq. (^0|)). Because of the commutation rule [Qi_i, Py] = 
iS^v, the CBH formula (|5^) makes the different fermions species anticommute. This must be properly 
extended to left-moving fermions, for instance by enumerating them consecutively after the right 
fermions (unless one uses the infinite- line quantization scheme). 

5.6 Bosonization of interactions 

Let us finally turn our attention to the basic interactions (|^). This requires introducing spin: we 
need two bosons Lp^ and Lpi . The boson fields may be naturally combined into spin and charge 
components: 

'Pc = —^{ip] ifs ^ —f={^] ~ fi) (100) 



and likewise for the chiral components 



Backward scattering 

The backward scattering term {gi) does not need normal ordering, since all four Fermi fields anti- 
commute. A straightforward application of the bosonization formulas is sufficient: 



Hi 



2 X! VaTlaTl-aV-a exp 



vgi 

(2^) 

vgi r _ _ 
(2^ h^T^^i 



(101) 



**In our periodic quantization scheme. On the infinite hne, with vanishing boundary conditions at a:: = ±oo, this is 
not true and [<ji(a;), ^(x')] = ^i. This garantees the anticommutation {'ip(x),'ii}{x')} = 0. 
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Since all the Klein factors present anticommute, '7|'7|?7x?7x = VlVlV^V] ^^'^ it remains 

Hi = ^V^mVlVl cos{V8^(ps) (102) 

The product of Klein factors poses a slight conceptual problem. The Klein Hilbert space (on which 
the 7]^ and fj^ act) is the representation space of a Clifford algebra, and has a minimal dimension 
which depends on the number of fermion species. With four species (two spins, times left and 
right), the minimal dimension is four, exactly like Dirac matrices in four-dimensional space-time. 
For instance, the following representation is possible, in terms of tensor products of Pauli matrices: 

?7| = CTi «) (Ti ?7| = 0-3 (g) fJi 

^1 = ® CTl ?7| = 1 (g) (72 



One readily checks that this is a faithful representation of the Clifford algebra (^8|). Then, the above 
Klein product is 

mv^vivi = i«)cr3 (104) 

This matrix is diagonal with eigenvalues ±1. If no other term in the Hamiltonian contains a nondi- 
agonal Klein product, one may safely pick a Klein eigenstate and stick with it, forgetting the rest 
of the Klein Hilbert space. This is some kind of gauge choice.^ In the case at hand, we may chose 
the eigenvalue -t-1 and write 

7^1 = i^cos(\/8^(^.) (105) 

Forward scattering 

The next interaction term (52) is quite simple to bosonize, since it involves only density operators: 

^ ^ d^ipadgipa' = '^^'^'^ d^ipcdzifc (106) 



and likewise for the spin interaction: 



n = ^^d,^,d,v, (107) 



Likewise, the forward terms (^4) are 



vg4:,c 

2w 



and 



Umklapp scattering 



[{d^iPcf + {d.^cf] (108) 



Finally, the Umklapp term (173) is treated in the same way as 171, this time with the charge boson 
involved: 



V93 

2(27r)2 



iv^vMi + vivmm) e^^^(*T+0i+*"T^-0i) + h., 

VI vim VI e'-^-^-^ +U.C. 



2V]VlV]Vlcos{v8T:ipc) (110) 



'One has to be careful when calculating various correlation functions, since the operators involved may also contain 
Klein factors, and simultaneous diagonality of all relevant operators with interactions is important if one wants to 
forget about the Klein Hilbert space. 
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After a mere transposition, the Klein prefactor is the same as (104). Therefore, 

^3 = ^ cos(V87r(y5c) 



(111) 



6 Exact solution of the Tomonaga-Luttinger model 

6.1 Field and velocity renormalization 

The power of bosonization hes mainly in its capacity of representing fermion interactions, namely 
7^2 and 7^4, as parts of a noninteracting boson Hamiltonian. The continuum electron model defined 
by the Hamiltonian density 



'Ht.L. = + 7^2, c + 7^2,5 + 7^4, c + 7^4,; 



(112) 



is called the Tomonaga-Luttinger model (TL model). The free electron part Tip may be expressed 
as the s um o f free boson theories for the spin and charge bosons, with the same velocity v. As seen 
in Sect. 5£, each interaction term is expressed either in terms of the spin boson or of the charge 
boson. Thus, the total Hamiltonian is the sum of charge and a spin Hamiltonians, i.e., the two 
sectors (charge and spin) are completely decoupled:^^ 



He + Ti-s 



(113) 



When expressed in terms of the spin and charge bosons ~ more precisely, in terms of the conjugate 
the Hamiltonian of each sector becomes 



field H^ and of dxffi 



Hfj. = 

^4,p — 



27r ^ 



\- H'^ where 

m + (9.<p^)'] 



where fi — c, s. Combining these expressions, one finds 



Hn 



1 



where 



K,. 



71" + 52,M + 54, M 



1 



.94,p 
TT 



(114) 



(115) 



(116) 



The constant is simply a renormalization of the field (pfj_. This is most obvious in the associated 
Lagrangian: 



L,. 



1 



2K„ 



dx 



1 



(117) 



The Hamiltonian (115) may be brought back in the canonical form ( pOD simply by introducing 
rescaled operators 

v'^ = -^v^, H' = ,/K;.n. 



(118) 



^"^This is true at the field level, but not completely true at the spectrum level, where things are more subtle, as 
shown above. Obviously the number of spin and charge excitations are not independent, since the total number of 
electrons is conserved. 
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which still obey the canonical commutation relations. Clearly, the field ip'^ has radius 

R^^^ = (119) 



Remarks: 



1. Let us point out that we have supposed the couplings 172. and f;4_^ to be momentum inde- 
pendent. However, this is not necessary: the above solution of the Tomonaga-Luttinger model 
may be repeated with momentum-dependent 32. /j (9) and 34.^1(9), for the simple reason that 
the boson Hamiltonian is quadratic and that different momenta are therefore decoupled. In 
particular, the above formula for the velocity renormalization may be used to write down the 
exact dispersion relation of charge excitations: 




^(fc) = = .|fcU/ ( 1 + ^liiM) - ( 3JhM\ (120) 



For simplicity, we shall ignore from now on such a momentum dependence of the couplings, 
mainly because it is generally irrelevant (in the RG sense). 



It is important that the perturbations 7^2 and H.^ of Eq. (114) be considered as perturbations 
of the Hamiltonian, not of the Lagrangian. If the conjugate momentum H were replaced by 
dt^p/v and the resulting expressions considered as perturbations of the Lagrangians, the field 
and velocity renormalization would be quite different. The reason is that, as they are expressed 
here, these perturbations depend on dt^p^ and would change the definition of the conjugate 
momentum, were they to be applied to the Langrangian. 



6.2 Left-right mixing 



The perturbations 'H2.11 of (114) mix right and left bosons. One then expects that the new eigenstates 
will be linear combinations of left and right excitations of the original boson. Indeed, rescaling the 
field ip> ((/?' = (p/^fK) implies the opposite rescaling of its conjugate momentum H - and therefore 
of the dual field i? - in order to preserve the canonical commutation rules: -d' = -d^/K (we drop 
the spin/charge index for the time being). Thus, the left and right bosons </> and (j) are not simply 
rescaled by if, they are mixed: 



i(^-^)^0' = i ("-^^-VXi?) (121) 



Expressing the old right and left bosons (0 and 0) in terms of the new ones (0' and 0') and defining 
^ as K = e^^, one finds 

= cosh ^ 0' -I- sinh ^ 0' 

= sinh ^ 0' cosh ^ 0' (122) 

In terms of the mode expansion (^o|), this left-right mixing is a Bogoliubov transformation of the 
creation and annihilation operators, if one collects like powers of e**^^ at equal times: 

bn — cosh^ h'^ -f sinh^ b'J 

In = sinh ^ cosh ^ 5; (123) 
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Thus, the fermion field ^/;|, a pure right-moving field in the unperturbed system, becomes a 
mixture of right- and left-moving fields in the Tomonaga-Luttinger model, if if^ ^ 1 (i.e., if gi ^^ ^ 0). 
Consider for instance the case Kg — Kc ^ 1: 

1 



V'T(a;,T) 



/2tt 
1 



7/1 e~'^2''"*'= e"''^-^ 



/27r 

iV27rcosh4 0^ g— iV27T- sinh { 0^ g— iV2iT-()is 



(124) 



This is expected, since the interaction 0,2 surrounds the free electron with a cloud of right and left 
electrons. 



6.3 Correlation functions 
Green function 



Once the renormalizations (116) are obtained, the calculation of physical quantities (correlation 

functions) is a straighforward task in principle, at least in real space. Let us start with the one- 
particle Green function 

GT(x,r) = (CT (a;, r)cj (0,0)) (125) 

Near zero wavevector, in terms of continuum fields, this becomesp] 

Gt ix,T) = (V-T (x, r)7A| (0,0)) + (i^T (x, t)^ (0, 0)) (126) 

Let us co ncentrate on the right-moving part, in the case Kg = 1 and Kc ^ 1 for simplicity. Using 
Eq. (124), it becomes 

(V'T(a;,T)V't(0,0)) = J_(c-'^^™^h«0^(2e)giv^ cosh? 0^(0)^ 

27r 

^ ^ giV27rsinh5 0^(2^) g— iV27rsinh? 0^(0)^ 
^^g-jV2^,/.',(z,) giv^^UO)^ 
11 1 1 



27r (w^T - z.x)^'=°''h=^« (ucT -f ia;)5s'"i^^C (u5r-ia;)i/2 
11 1 1 

27r (wcT — ix)^!"^ \vrT — ia;|*'= (wsT — zx)^/^ 



(127) 



where the velocity difference between spin and charge demands a distinction between spin and charge 
complex coordinates {zg and z^). We have introduced the exponentp] 

^c- J(i^c + ;^-2) (128) 

Once the Fourier transform is taken in time and space, this Green function leads to extended spectral 
weight,^ as illustrated on Fig. |[ This spectral weight has power-law singularities near uj = ±Vcq 
and Lo — VsQ, with exponents related to Oc- From this spectral function, it can be shown that the 
momentum distribution function n{k) has an algebraic singularity at the Fermi level: 

n{k) — n{kF) — const, x sgn(fc — fci?)|fc — kpf" (129) 

and so does the one-particle density of states: N{uo) ^ \^\^''- 

^^There is also a singular component near k = 2kp is ii" < 1. 
-"^■^The notation 7c = \dc is also used, as is = dc- 

^^For a detailed calculation of the spectral weight from the above space-time Green function, see Ref. [p| 
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Figure 5: Typical spectral function A{k,Lu) of the Tomonaga-Luttinger model, with Ks = 1. The 
exponents appearing on top refer to the three singularities: {lj + Vcq)^", {lj — Vsq)^"^^^^ and {uj — 



Compressibility 

Next, let us consider the compressibility near q — 0: 

x{x,t) = {ntot.{x,T)ntot.{0,0)) 

Since n^^tot. = dxifi^/y/n, we have 

[2 [T 
"tot. = \/ -dx^fc whereas 5^ = \ —dx(ps 

V TT \ ZTT 

is the z component of the uniform spin density. The compressibility is then 

2 



(130) 



(131) 



-{dx(Pcix,T)dx(PciO,0)) 
TT 

9K 

TT 

IK 

= [{d.v'Ax, r)d.v^'M 0)) + {d-.v'M r)d,yM 0))] 

TT 

Kc_ 
27r2 



(132) 



This, of course, is the quasi-uniform component of the compressibility. The corresponding result for 
the magnetic susceptibility (still at A% = 1) is 



(5.(X,T)5.(0,0)) = ^ 



1 



(133) 



where Vg is the spin velocity. 

The compressibility also has a charge-density-wave component at ±2fci?, determined by the cor- 
relations of the operator 



CDW 
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K,-\ 



Figure 6: RG flow with backward scattering gi near the point Kg = 1. 



27r 
1 

= Y' 

The CDW compressibility is then 
(OcDw(a:,r)O^Dw(0,0)) = 



-2ikFX - -i\/iT!ipa 



(134) 



27r2 



g -2ikpx ^ ^iV2TTipc{z) g-iV27r(^c (0) ^ ^ giV27rips (z) g-iV27r(^s (0) 



-2ikFX 



27r2 



l^cT — ia;|-^'= jusT — ix\ 



(135) 



Notice that the Klein factors disappear from the correlation function, as they do in general from a 
product of the type OOK 



6.4 Spin or charge gap 

The interactions TLi and Ti,^ do not lend themselves to an exact solution like 7^2 and 7^4 do. However, 
they are subject to RG arguments, at least at weak coupling. 

Let us first consider the backward scattering term (105). Combined with the spin boson Hamil- 
tonian (for Ks ^ 1), the perturbation defines the well-known sine-Gordon model, one of the most 
studied one-dimensional field theoriesiR 



1 



UsG = [Uf + (a,^'j2] -t- ^ cos ( V&^sV's 



27r2 



(136) 



(we used the rescaled boson ip^ in order to give the unperturbed part its standard normalization). 
If Ks = 1, the cosine term is marginal, as can be checked from Eq. (^7|): A = h + h = l + l = 2. 
One would naively conclude that it is irrelevant if Kg > 1, and relevant if i^s < 1, thus leading to a 
gap in the spin excitations. The story is a little more complicated by the existence at 51 = not of 
a fixed point, but of a fixed line, parametrized by Kg. It turns out that Kg is subjected to RG flow 
away from that line and one must instead consider the set of coupled RG equations 



d£ 



1 



= -2gi{Kg - 1) 



(137) 



first obtained by Kosterlitz and Thouless jg]. These equations are demonstrated in Appendix A, 
from the short distance product, or OPE, of the perturbations involved. The RG flow is illustrated 



^*For the sake of argument, I have neglected the velocity renormalization here {vs = v). 
additional renormalization of gi, defined by vgi = Vsg[, which is neglected at low coupling. 



It would lead to an 
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on Fig. ^. One checks that, at weak coupling, the system flows towards strong couphng if \gi] > 
2tt(Ks — 1) is then said to be marginally relevant) while it flows to gi = if \gi] < 2tt(Ks — 1) 
(gi is then said to be marginally irrelevant). Thus, we expect a spin gap for gi infinitesimal and 

This expectation is fuUfiUed exactly at the special value Ks = \. There, the cosine term has 
scaling dimension 1 and is precisely what one would obtain from bosonizing the Dirac mass term of 
Eq. (|2^), except that the corresponding fermion represents not electrons, but massive spinous. This 
Ks = 1 theory is called the Luther-Emery model Q and is a free theory for spinous, even though it 
is a complicated interacting theory for electrons (and who knows what happens in the charge sector 
at the same time. . . ). 

The above analysis can be repeated for the Umklapp term (111), which occurs only at half filling. 
One simply has to replace gi by g^ and Ks by Kc- Thus, a charge gap develops if K,, < 1 (i.e. if 
52, c is repulsive) or if jgs] > 2tt{Kc — 1). Again, free massive holons occur at Kc — \ - 



7 Non-Abelian bosonization 

7.1 Symmetry currents 

The bosonization procedure described in Sect. ^ is often called Abelian bosonization, because a 
single compactified boson cp has a C/(l) symmetry if ip + a and the group U{1) is Abelian (i.e. 
commutative). Abelian bosonization may be used to describe the spin sector of a ID electron gas, 
but this description is not manifestly spin-rotation invariant. Indeed, the chiral currents associated 
to the three spin components are easily shown to be 

J' = ^V4(c^3)a/3V'/3 

J"" = ^V4(c^l)a/3V'/3 
= ^i^i{<^2)a(3i^P 

The z component manifestly has a special role in this representation. While the correlations 
{J°'{x)J'^{0)) (a — x,y,z) all decay with the same power law, they do not have the same nor- 
malization, a somewhat unsatisfactory feature. 

Fortunately, there exists a bosonization procedure that is manifestly symmetric under the spin ro- 
tation group SU(2), or under any Lie group for that matter. This non- Abelian bosonization expresses 
a set of fermion fields in terms of a matrix field belonging to a representation of a Lie algebra, in- 
stead of one or more simple boson fields. The corresponding boson theory, the Wess-Zumino-Witten 
(WZW) model, is far less familiar than the ordinary boson. But in all situations where symme- 
try considerations (including symmetry breaking) are important, non- Abelian bosonization is the 
method of choice. It is of course unrealistic to give a rigorous and complete introduction to the 
WZW model and non- Abelian bosonization within this brief tutorial. Ref. ||] may be consulted for 
an in-depth discussion. 

Let us first return to Abelian bosonization and insist on the role of the U(l) current, before 
generalizing to larger symmetry groups. The currents associated with the U(l) symmetry are given 
in Eq. ( ^o[ ) for a single species of fcrmions and coincide with the chiral fermion density fluctuations. 
For the charge U(l) symmetry, it is the total (t+i) density fluctuation that matters, and the charge 

^^Note that the sign of gi is of no consequence: it is the result of a gauge choice through Klein factors. The 
arbitrariness in the sign of gi may be traced back to its expression in terms of left and right fermions: changing the 
sign of either or i/ij is of no consequence in the noninteracting case, but changes the sign of gi and 33, 

without affecting those of 32 or 34 . 



dzips 



2tt 



mVi sin (^\/87r0< 



-^V'iVl cos ( V»nq 



(138) 
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currents are therefore 

ii — % 

Jc = —f^(^z'Pc Jc = -=dz(pc (139) 

v27r v/Stt 



In terms of the U{1) group element g — e^^' , these currents are proportional to g ^dzQ and g ^dgg. 
The charge Hamiltonian can then be expressed in terms of these currents, as in Eq. (B4): 



H^^Jdx {Jl + Jl) (140) 

This is enough to define the theory, provided that (i) the commutation relations of the currents are 
specified: 

[,h{x),J,{x')]=~-5'{x-x') (141) 

TT 

and (ii) the vertex operators are introduced, with their proper commutations with the currents. 
Note: The above current commutation relation can be demonstrated in the same way as we have 
recovered the fermion anticommutation relations (^7|) , from the current-current correlation function 

obtained from Eq. (|59| ) by differentiating w.r.t. z and z' . 

Likewise, non-Abelian bosonization (for SU(2)) may be introduced starting with the various 
components of the spin currents J a defined above in terms of two electron fields V'o-, and expressed 
in terms of a group element g as 

J{z) = {d.9)g-^ = \Y. •^-'^^ •^(^) = 5"'(5.-3) ^\Y. JaCJa (143) 

a a 

The Hamiltonian of the SU(2) WZW model can be expressed in terms of these SU(2) currents: 



o 3 



^ y da; [.r(x)r{x): + : /"(x) J"(x) :] (144) 



where the integer k is called the level of the WZW model, and where the currents now obey the 
following commutation rules 

[r{x),j\x')] = -~'^5ab6'{x-x')+ieabeJ'{x)5{x-x') 
ZTT 

[r{x)J\x')] - -'A5ab5'{x-x')+ieaber{x)5{x-x') 
ZTT 

[r{x),j\x')] = (145) 
Once the currents are Fourier expanded as 

r[z) ^^^Jn e'2™^/^' (146) 

n 

The above commutators translate into the following set of commutation rules, known as a Kac- Moody 



ra: 



[J^Jn] = iSabeJ^+n + ^kndabSm+n,0 (147) 
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States of the WZW model are created by applying the operators J° {n < 0) on the vacuum |0) and 
on a finite number of spin states. 

The formulation of the WZW in terms of currents is not as fundamental as its definition in the 
Lagrangian formalism, in terms of a field g{x, t) taking its values in SU(2) (or in another Lie group). 
The model must be such that the currents are chiral, i.e. dzJ — dgJ — 0, and must have conformal 
invariance. The needed action is 

S[9] = 1^ /d'^ trid'^g-'d.g) j d^x e^^xtv {g^' 8^ g g-' d'' g g-' d"" g) (148) 



B 



where the second term is topological: it is integrated in a three-dimensional manifold B whose 
boundary is the two-dimensional space-time, but its value depends only on the field configuration on 
the boundary of B. The level k must be an integer in order for the WZW model to have conformal 
invariance and chiral currents. Working with the action ( |14^ ) is somewhat unwieldy and almost 
never done in practice, once key identities (Ward identities) have been derived, leading, among 



others, to the commutators (145) 



Bosonization cannot be done in terms of the currents alone. The group elements g occur in the 
representation of left-right products of fcrmions, as specified by the Witten formula^^ : 

' Zn 

Ml = T^(3^)o/5e^^^= (149) 
Zn 

where the presence of the charge boson tpc is still necessary in order to represent the charge component 
of the products. This formula leads to the following representation of the charge and spin densities: 



27rntot. (x) = J + J - 



e-2»'=^^ e^^-^Hrg -I- H.c. 



2ttS^^^{x) = J" + J''+ [e-2'fc^-^e*^'^<^tr((7(Ta) + H.c.] (150) 

The WZW model is a conformal field theory, with a set of well-identified primary (or scaling) 
fields. It is completely soluble, in the sense that the structure of its Hilbert space is completely 
known and correlation functions obey linear differential equations. The central charge of the level-fc 
SU(2) WZW model is 

'=^2 ^'''^ 
and it contains scaling fields of all spins from s = to s = A:/2, with conformal dimensions 

k + 2 ^ ^ 

The spin part of the ID electron gas with one band is described by the simplest of all WZW models: 
SU(2) at level 1. The matrix field g corresponds to s = i and has conformal dimensions (i, j). Its 
components may be expressed in terms of the spin boson ips, if we follow Witten's formula and the 
Abelian bosonization formulas: 

The equivalence of the k — 1 SU(2) WZW model with the spin boson theory is possible because the 
k = 1 model has central charge unity. Such a correspondence between a WZW model and simple 
bosons is possible only if c is an integer. 
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Other WZW models describe more exotic critical systems. The biquadratic spin-1 chain, with 
Hamiltonian 

if - ^ S„ • S„+i + /3(S„ • S„+i)2 (154) 

n 

is critical at /3 = 1 and /? = — 1. At (3 = —1, the low-energy limit is described by the k — 2 SU(2) 
WZW model, which has central charge c = |. This model may be expressed in terms of three 



Majorana fermions |11, Ij. At /? = +1, the low-energy limit is described by the fc = 1 SU(3) WZW 
model. A recent attempt at extending Zhang's SO (5) theory of antiferromagnetism and d-wave 
superconductivity to one-dimensional systems is formulated in terms of WZW models |0 . 



7.2 Application to the perturbed Tomonaga-Luttinger model 

The noninteracting part of the Tomonaga-Luttinger Hamiltonian may be expressed as follows in 
terms of charge and spin currents: 

Ho = f [J.^ + J^]+^[J^+P] (155) 

It is the sum of a free charge boson and of a A: = 1 SU(2) WZW model. The various interactions 
bosonized in Sect. 5.6 have the following expressions in terms of currents: 

Hi = -2vgi{JxJx + JyJy) 
HI = \vgi,c[Jl + Jl] 

^2 = Avg2,sJzJz 

Hs = ^cos(V8^(pc) (156) 

This correspondence may be established with the help of Eq. ( |138[ ). It is then manifest that 34 c 
and gi^s only renormalize the velocity. Moreover, if gi = — 2(72,s, and 7i| combine into a single, 
rotation-invariant interaction: 

Hi+n^^ ivg2,s3 - J (51 = -2g2,s) (157) 



The spin sector 

Let us restrict ourselves to this isotropic case, at half-filling, so that 173 ^ 0. As we have seen in 
Fig. ^, the RG flow will bring g^ to strong coupling if jgaj > —2g2^c and the charge sector will have 
a gap. At the same time, the coupling gi flows along the separatrix of Fig. ^ towards the fixed 
point gi — 0, Kg = 1 (a detailed RG analysis, following the method described in Appendix A, shows 
that the separatrix is precisely determined by the condition of rotation invariance gi — —2g2.s)- 
Thus, the low-energy limit of the half-flUed perturbed TL model is the k = \ WZW model, provided 
I53I > -252, c and gi = -2g2,s- 

The appearance of a gap in the charge sector has another consequence, that of changing the 
scaling dimensions of some operators. Specifically, let us consider the 2kp component of the spin 
density, which becomes the staggered magnetization at half-filling: 

ST^ = iV^l(a,),,^,+H.c. 

- -^5..K).,e^^^=+H.c. 

= --^tr(aaff)e'^'^=+H.c. (158) 
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where we have appHed Witten's formula (|149|). The cosine Umklapp term gives Lpc a nonzero 
expectation value, such that (cos ^/AiTLpc) ^ and (sin ^/^Lpc) = 0. Neglecting fluctuations around 
that expectation value, one may write the staggered magnetization density as 

Sl^" ^ const. (tr(CT„5) + H.c.) (159) 



According to Eq. (152) with s = ^, this operator has scaling dimension A = /i + /i = i, although it 
had A = 1 initially (before the onset of the charge gap) . Essentially, the charge contribution to the 
fluctuations, i.e. the operator e**^'^'^=, also of scaling dimension A = i, is frozen; this makes the 
operator S'^^'' more relevant. The staggered magnetization correlation is therefore decreasing like 
1/r, instead of the 1/r^ decay of uniform correlations: 

(5f-(x,i)5f-(0,0)) ^ (160) 

Other Perturbations 



One may add various relevant or marginally relevant perturbations to the model (155): 



1. An Ising-like anisotropy would contribute a term of the form Tiising = ^JzJz- This term, a 
correction to 172,5 is marginally relevant if A > 0: it brings the flow above the separatrix of 
Fig. ^. On the other hand, it is marginally irrelevant if A < 0, and brings the system towards 
an XY fixed point. 

2. An explicit dimerization, maybe in the form of a staggered hopping term in the underlying 
Hubbard model. This perturbation would be proportional to the CDW (i.e. 2kF) component 
of the electron density. In terms of the WZW field g, this may be obtained by replacing the 
Pauli matrices of Eq. ( |158D by the unit matrix. At half-filling, the freezing of cos ^/Airipc would 
leave the single scaling field tvg in the perturbation: 

Hdim. - Atr5 (161) 

This term has scaling dimension A = ^ and would lead, according to Eq. (p^), to a spin gap 
of order m ^ A^/'^. 

3. A uniform external magnetic field B. This couples to the uniform magnetization density Jz+Jz, 
producing a perturbation with nonzero conformal spin. It breaks rotation invariance and is 
best treated within Abelian bosonization. In terms of the spin boson cps, the perturbation is 
Xdx^s, where X oc B. The spin Hamiltonian is then 

n^^v [nl + (d^ips)^] + xd^^s (162) 

(we have supposed Kg — 1, i.e., rotation invariance in the absence of external field). This 
perturbation can be eliminated from the Hamiltonian by redefining 

ips ^ '■Ps- -X (163) 

V 

In terms of the original electrons, this is almost equivalent to changing the value of kp (almost, 
because it should have no effect on the charge sector). Basically, the uniform magnetic field 
introduces some incommensurablility in the spin part of the system. 
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8 Other appUcations of bosonization 

8.1 The spin-| Heisenberg chain 

The spin-i Heisenberg model in one dimension may also be treated by bosonization in the low-energy 
limit. The microscopic Hamiltonian is 

H = ^ {j{Sf^S^^i + S^^S^^i) + JzS^Sji^i} (164) 

n 

This Hamiltonian may be mapped to a spinless fermion problem via the Jordan- Wigner transfor- 
mation: 

= 4 CXp ^i-K 4nCm^ = C^^Cn - ^ (165) 

where the operators c„ and cj„ obey standard anticommutation relations. In terms of these operators, 
the Heisenberg Hamiltonian is 



-H" = XI + H.c.) + Jz{cic„ - ^)(cl+ic„+i - i)| 



(166) 



The fermion number cJjC„ is related to the total magnetization M, in such a way that half-filling 
corresponds to M — 0, the appropriate condition in the absence of external magnetic field. The 
continuum limit may then be taken and the bosonization procedure applied. The spinless fermion 
is then endowed with interactions proportional to Jz . Two of these interaction terms are of the (72 
and (74 type and renormalize the boson field ip and the velocity v. From comparing with known 
exact exponents, one derives the following relation between the radius R of the boson and the bare 
interaction strength J^: 

i?=-^(l-icos-i f^H (167) 



Various correlation functions may be calculated, once the correspondence between spin components 
and the boson field is known : 

S'ix) = -l-a,^H-(-l)^cst.e'^/« 
ZttH 

For instance. 



1 

16"^ i? 

+cst.(-l)^(e*'^(="'*)/-«e'"^(°'")/-") 



{S'{x,t)S' (0,0)) = ^^—^{{dz^{x,t)dzV^{0,0)) + {dMx,t)d-MO,m 



1 



1 1 



(x - wi)2 (a; + vif 



Note that the staggered part has a nonuniversal exponent, which depends on Jz/ J ■ At the isotropic 
point Jz — J, R = l/-\/27r and this exponent is ^, exactly like the non-Abelian bosonization 
prediction (160). 
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Figure 7: Fluctuation of the edge of a quantum Hall system. 



8.2 Edge States in Quantum Hall Systems 

The two-dimensional electron gas (2DEG) in a magnetic field has been the object of intense theo- 
retical investigation since the discovery of the quantum Hall effect (integral and fractional) in the 
early 1980's. While the integer Hall effect may be understood (in the bulk) on the basis of weakly 
correlated electrons, the fractional quantum Hall effect is a strongly correlated problem. Quantum 
Hall plateaus correspond to gapped states, but electric conduction nevertheless occurs through elec- 
trons at the edge of the 2D gas, the so-called edge states. One must imagine a 2DEG immersed in a 
magnetic field B and further confined to a finite area by an in-plane electic field E (in, say, a radial 
configuration). The field causes a persistent current along the edge of the gas: 

j = a.^yi A E (170) 

where Uxy is the Hall conductance, equal to ve^, where v is the filling fraction. The edge electrons 
drift in one direction at the velocity v = cE/B. We will adopt a hydrodynamic description of the 
edge, following Ref. |0. Let x is the coordinate along the edge and h(x) the transverse displacement 
of the edge with respect to the ground state shape of the gas (cf. Fig. If n = v/2ttI\ is the 
electron density [is is the magnetic length), then the linear electron density in an edge excitation 
is J{x) — nh{x). The energy of the excitation is the electrostatic energy associated with the edge 
deformation: 

H = 1 dx ^e^hJE = TTiyv I dx j2 (171) 



The field J being from the start treated as a collective mode, bosonization is natural. The main 
difference with the ID electron gas is that only one half (left of right) of the boson is necessary, since 
the edge excitations only propagate clockwise or counterclockwise, depending on the field direction. 
We therefore assert 

j = ^a,0 = -^a,0 (172) 

and the Hamiltonian density for edge excitations is 

n = -vv{d,(j)f (173) 

The filling fraction v plays the role of a field renormalization K = 1/v. The edge electron must be 
represented by a field ^j{x) such that 

[J(x), V(a;')] = -5{x - x')i}{x) (174) 
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From the short-distance product of J with e one finds that the correct representation for V' is 

^(x) = e^VI^"^/- (175) 
V 27r 



From this we derive the one-electron Green function 

Note that this coincides with a free-electron propagator in the case v = 1 only, i.e., when the first 
Landau level is exactly filled. Even in the integer effect, the edge system is a Luttinger liquid if 



TO = 1/v > 1, albeit a chiral Luttinger liquid. In momentum space, the Green function (176) is 



Lo ^ vk 

and the momentum-integrated density of states near the Fermi level is 

iV(w) (X l^r^i (178) 

8.3 And more. . . 

This short tutorial has concentrated on basic Luttinger liquid physics. Time and space are lacking 
to introduce the many physical systems that have been treated with bosonization. Here we will just 
indicate a few useful references. 

The bosonization formalism may also be used to treat disorder. The problem of an isolated 
impurity in a Luttinger liquid is treated in Refs. |^ and |l^, following Kane & Fisher's original 
work JlGf . The more complicated problem of Anderson localization, resulting from the coherent 
scattering off many impurities, is also described summarily in Ref. JTsf . Formally related to impurity 
scattering is the question of electron-phonon interaction in a Luttinger liquid, a review of which may 
be found in Ref. ||]. 

The Kondo problem, or the strong interaction of magnetic impurities with conduction electrons, 
has also been treated using bosonization, and a strong dose of conformal field theory, by Affleck and 
Ludwig. Ref. ijl^ is a comprehensive review of the progress accomplished with the help of conformal 
methods. Most interesting in this application is the three-dimensional nature of the initial problem, 
before it is mapped to a one dimension. 

Spin chains of spin s > ^ can be constructed from 2s fermion species with a Hund-type coupling. 
This is the approach followed in Ref. also reviewed in Ref. [Q. The spin-1 chain has been 

treated as a perturbed fc = 2 SU(2) WZW model in Ref. 

Coupled spin chains, or spin ladders, have been studied in Abelian bosonization and are reviewed 
m Ref. 131. A non- Abelian bosonization study of two coupled spin chains is performed in Ref. W. 
Doped spin ladders are coupled Luttinger liquids. Bosonization studies may be found in Refs. |]l9| 
and ll^. Ref. jl^ looks at this problem from the perspective of S0(5) symmetry. The latter two 
references treat the question of d-wave fluctuations in a doped ladder. 

The reader is referred to Ref. ||l^, js) for additional references (we apologize for our poor bibli- 
ography). More applications of bosonization to condensed matter can be found in Ref. Other 
important review articles on the subject include Refs. and A historical perspective is given 
in the collection of preprints ^ . 

9 Conclusion 

Let us conclude by highlighting the principal merits and limitations of bosonization: 
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1. Bosonization is a nonperturbative method. The Tomonaga-Luttinger (TL) model, a contimimn 
theory of interacting fermions, can be translated into a theory of noninteracting bosons and 
solved exactly. 

2. Bosonization is a method for translating a fermionic theory into a bosonic theory, and even- 
tually retranslating part of the latter into a new fermionic language (cf. the Luther-Emery 
model or the use of Majorana fermions). This translation process is exact in the continuum 
limit, but does not warrant an exact solution of the model, except in a few exceptional cases 
(e.g. the TL model). For the rest, one must rely on renormalization-group analyses, which 
generally complement bosonization. 

3. The bosonized theory may have decoupled components, like spin and charge, which are not 
manifest in fermionic language. Thus, spin-charge separation is an exact prediction of bosoniza- 
tion (even beyond the TL model). When one of the components (spin or charge) becomes 
massive, this favors the emergence of physical operators with new scaling dimensions, like the 
staggered magnetization of Eq. ( |159| ) . 

4. The contact with microscopic model is not obvious, since an infinite number of irrelevant 
parameters stand in the way. But bosonization leads to universal or quasi- universal predictions, 
and one can argue that such predictions are preferable to the exact solution of a particular 
microscopic Hamiltonian of uncertain relevance. 

5. In the rare cases where a microscopic Hamiltonian has an exact solution (e.g. the ID Hubbard 
model), the solution is so complex that dynamical quantities cannot be explicitly calculated. 
Bosonization can then be used as a complementary approach. Exact values of thermodynamic 
quantities of the microscopic theory may be used to fix the parameters (e.g. Uc, Ws, Kc, Kg) 
of the boson field theory, and the latter may be used to calculate dynamic quantities (e.g. the 



spectral weight) or asymptotics of correlation functions p^, 22 



6. The low-energy limit corresponds to a few regions in momentum-space where bosonization has 
something to say: the vicinity of wavevectors that are small multiples of kp. 

7. Finally, bosonization is limited to one-dimensional systems, despite attempts at generalizing 
it to higher dimensions H, E3| . 
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A RG flow and Operator Product Expansion 

In this appendix we show how to derive the approximate RG trajectory around a fixed-point from the 
Operator Product Expansion (OPE) of the marginal perturbations involved. Consider a fixed-point 
action Sq describing a conformal field theory, and a set of perturbations with coupling constants gi: 

S[^] = Som + J2g^ f dxdTO,{x,T) (1^9) 

Here $ denotes a field or a collection of fields that defines the theory, in the sense of the action 
and of the path integral measure. The local operators Oi are in principle functions of The RG 
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procedure may be performed in momentum space or in direct space. It is more convenient here to 
used the direct space approach and to divide space-time into tiles of side L. A space-time point may 
then be parametrized as 

x = X + y (180) 

where x is a space-time point, X is the space-time coordinate of a tile and y a coordinate within 
a tile. The RG step is performed by tracing over degrees of freedom labeled by y, leaving us with 
"block-spin" degrees of freedom $(X). This is done approximately be expanding the perturbations 
to second order in the partition function, in order to construct an effective action for the "block 
spin" variables: 

exp-5cff.[$] - e-^«[*l |l - Eff'^ / dx(0.(x)>> + ^^5.5^/ dxdx'(0,(x)0,(x'))> + • 

(181) 

where the expectations values are taken over the fast modes, i.e., the modes that vary within a tile, 
in the fixed-point theory. Notice that the fixed-point action 5*0 is the same for $ as for the block-spin 
field $, precisely because 5*0 represents a fixed point. 
Consider first the integral 

/'dx(0,(x))> =E / dy(0,(X + y))> (182) 
X ^^^^ 

The integration over the tile of side L may be inferred from dimensional analysis: if {hi, hi) are the 
conformal dimensions of the operator Oi and Ai = hi + hi is the corresponding scaling dimension]^ 
then this integral must be (there is no other quantity with the correct scaling dimension, 

since the correlation length of the fixed point theory is infinite), times some effective X dependence, 
so that the first term of Eq. (|l8l|) may be written as 



i X 

where the field Oi(X) denotes this X dependence. 



The second term of Eq. (181) may be treated in the same way, once the product of operators has 



been reduced to a sum of single operators by the OPE: 

0.(x)0,(x') -E^I^O^("') |^,x.|i.+A,-A. (184) 

It is assumed here that the set of operators Oi forms a "closed algebra" under this short-distance 
expansion. Only the most divergent terms really matter. Let us substitute this OPE in the second 



term of Eq. (181) and evaluate the integrated expectation value like we did for the first term: we 
find L^^^*-', times 0(X), times the integral over the relative coordinate x — x', which scales like 
2^2-A,-A,+Afc^ or fike InL if 2 = Ai -I- Aj - A^. The effective action, after this RG step, is thus 
given by 

exp-^cff.i-l] - e-^«[*l|l-E5fcydXL2-A.o^(X) 



1 r r2-Ai-A,+Afc ~l f 

\T.T.s^9,^^ }^^-^7dxo.(x) 

fc i,i 



(185) 



^^For simplicity, we will deal with operators of zero conformal spin only, i.e., we suppose that hi = hi. Operators 
with nonzero conformal spin naively do not contribute: they have zero expectation value. 
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(a sum over tiles is equivalent to an integral over X, since L is the new 'lattice spacing'). After 
reinstalling the perturbations in the exponential, one may define a renormalized coupling g'f, as 



2-A,, 



i,3 



InL 



L 



2-Afc 



(186) 



This RG step transformation must be made infinitesimal in order to be translated into an RG 
flow equation. To this end, we set L very close to unity, in units of lattice spacing, a somewhat 
formal procedure. In the case of a nonmarginal interaction (A/j 7^ 2) the first term dominates in a 
perturbative sense and, after setting L — with £ small, we find 



5^ « fffc [1 + (2 - A,)^] ^ ^^(2-A,)5fc 



(187) 



as expected. A more interesting case is that of a set of marginal interactions (A^ = 2), where 

9'k=9k-lY.(^%9^9,l^-^L ^ ^^AY^c%g^g^ (188) 

id 



2 ^ --v^ 



Let us apply this last equation to the spin sector of the ID electron gas subjected to the interac- 
tions gi and 52, in order to recover the flow equations (137). Let us first note that, when expanded 
to first order in 32. s from Eq. (|116|), those fiow equations become 



-9192 

TT 



dg2 
d£ 



1 2 
Zn 



(189) 



(we dropped the 's' index from now on: it is understood that we deal with the spin sector). The 
corresponding operators are 



Oi 



1 



cos(v 87ri^) 



O2 — —d(pd(p 

TT 



(190) 



The corresponding OPEs are obtained by applying Wick's theorem. Their most singular terms are 



Oi{z,z)Oi{w,w) 

02(2, z)02{w, w) 
02(2:, z)Oi{w,w) 



1 



1 



1 



167r4 (z-w)2(z-u;)2 27r3 

— \ 7 ^ rr dipBip 

TT'* [Z — W){Z — W) 

1 1 1 



1 



{z — w)' 



1 



{z — w)-^ 



1 



AtT^ (z — w)"^ {z — w)'^ TT^ \_{z — w)' 
1 1 r / N 

cos(v 87riy9) 



{z — w)' 



27r4 (z — w){z — w) 
This last OPE is derived from the relation 

9(p(z) e*"^("'''^) 



(191) 



ia 1 
At: z — w 



(192) 



which is obtained by applying Wick's theorem to the series expansion of e^°"^ (see, e.g., Ref. 
We may rewrite these OPEs as 



Oi{z, z)Oi{'w, w) 

02(2, z)02iw, w) 

02(2, z)Oi{w, w) 



const. 



1 



1 



27r-^ (z — w){z — w) 
const. + conf. spin 
2 1 



02(10, w) + conf. spin 



TT^ (z — w){z — w) 



Oi{'w,w) 



(193) 
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where 'conf. spin' denotes terms with conformal spin, which do not affect the beta functions as 
calculated above. From Eq. ( |188| ), we may then write the flow equations : 

dgi 1 1 2 nnA\ 

-di ^ -dl ^ 4^^i 



These differ from Eq. (18£) only by a rescaling of the flow variable £ by 27r, caused by our crude 



treatment of kinematic integrals. A corrected version of Eq. (18S) would then be 



-nY^C^g^g, (195) 



^,3 
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